OPEN 3 ACCESS Freely available online 



•0-PLOS I ONE 



Conservation Priorities when Species Interact: The 
Noah's Ark Metaphor Revisited 



Pierre Courtois, Charles Figuieres*, Chloe Mulier 

INRA-LAMETA, Montpellier, France 



Abstract 

This note incorporates ecological interactions into the Noah's Ark problem. In doing so, we arrive at a general model for 
ranking in situ conservation projects accounting for species interrelations and provide an operational cost-effectiveness 
method for the selection of best preserving diversity projects under a limited budget constraint. 

Citation: Courtois P, Figuieres C, Mulier C (2014) Conservation Priorities when Species Interact: The Noah's Ark Metaphor Revisited. PLoS ONE 9(9): e106073. 
doi:10.1 371/journal.pone.01 06073 

Editor: David L. Roberts, University of Kent, United Kingdom 

Received April 16, 2014; Accepted July 28, 2014; Published September 2, 2014 

Copyright: © 2014 Courtois et al. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits 
unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited. 

Data Availability: The authors confirm that all data underlying the findings are fully available without restriction. All relevant data are within the paper and its 
Supporting Information files. 

Funding: This study was supported by the ONEMA program 201 3-201 5. The funders had no role in study design, data collection and analysis, decision to publish, 
or preparation of the manuscript. 

Competing Interests: The authors have declared that no competing interests exist. 
* Email: figuiere@supagro.inra.fr 



CrossMark 



Introduction 

Weitzman [13] is a milestone in the economic theory of 
biodiversity. His "Noah's Ark Problem" is not only a modeled 
metaphor that is helpful to organize thinking on how to face 
conservation trade-offs with finite resources. It also results in a 
practical cost-effectiveness methodology that can serve as inspira- 
tion to guide conservation policies. The idea is, for each species i, 
to collect information about: i) C,, the cost of its protection, ii) 
APj, the increase of survival probability resulting from it, Hi) [/,, 
the direct utility of how much we value the species, iv) Dj, its 
distinctiveness. From this information, each species is assigned a 
number R' via the formula: 

AP 

Bf—^iPt+Ui), (1) 

which indicates its rank in conservation priorities. This ranking 
criterion has a theoretical foundation: it is rooted in a rigorous 
optimization model ([13], Theorem 4, p. 1295). 

This criterion sheds light on real biodiversity issues and has 
actually been used in several applications. Some of these have led 
to changes in allocation of conservation funding (e.g., in New 
Zealand; [9]), and variants have been used to allocate surveillance 
effort over space (e.g., [8]). Other applications are quoted in [5], 
But it is fair to say that this approach is more appropriate for ex 
situ conservation projects - say to build a gene bank or a zoo - 
rather than to manage a set of interacting species in their natural 
habitats. This is so because formula (1) uses no information of any 
kind about the web of life. Yet, in ecosystems, species interact. 
Some of them compete to share common resources, others develop 
synergies and mutually enhance each other or they simply pertain 
to the same trophic chain. Suppose, then, that the conservation 



authority has information about those ecological interactions, even 
if it is only under the rudimentary form of survival probability 
interdependencies. That is, it knows that a marginal increase of 
survival probability of species j will have an impact r,y on the 
survival probability of species i. Could this information be used to 
qualify formula (1) and increase its relevance when it comes to in 
situ conservation trade-offs? 

To our knowledge, three recent articles stress the need to 
account for ecological interactions in Weitzman's diversity 
concept. They have in common: i) to take into account the 
ecological interactions via interdependent survival probabilities in 
a simplified version of the Noah's Ark metaphor with two species 
[1], [11] or three species [12], ii) to show that this consideration 
can reverse the conservation priorities. The key of this note is to 
provide a general analysis of in situ conservation problems 
considering interdependent survival probabilities. Revisiting 
Weitzman's optimization problem, we extend his model in order 
to incorporate species interactions. Our principal output is to 
forward a general ranking formula that could be used as a rule of 
thumb for deciding in situ conservation priorities under a limited 
budget constraint. 

The sketch of the paper is the following. Section 2 incorporates 
ecological interactions in Weitzman's parable of Noah's Ark, with 
any arbitrary number of species. The crux of the section is to 
provide with a new rule for establishing in situ conservation 
priorities through the expression (12) below that encompasses 
formula (1) as a special case. The link between this formula and 
Noah's optimal policy is explained. Section 3 illustrates the 
relevance of this new formula within a two-species example. We 
check the robustness of our formula and end the paper with a 
discussion on the possibility of ranking reversal in relation to three 
stylized kinds of ecological interactions: predation, mutualism and 
competition. 
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Analysis 

The "Noah's Ark Problem" is a parable intended to be a kind of 
canonical form representing how best to preserve biodiversity 
under a limited budget constraint. In the initial version of 
Weitzman's modeled allegory, Noah's decision problem is, for 
each species i, to choose a survival probability between a lower 
and an upper bound, P, e [P,-,P/] , in order to maximize the sum 
of the expected diversity function: 

w({Pi}U), 

and the expected utility of the set of species: 



tf({p,}? = ,) = £ u >* p >- 
i=i 

Weitzman devotes much of his paper to defining the expected 
diversity function w({P;}f = i^ and to explaining its link with the 

concept of information content (see his Theorem 1, p. 1284). This 
function could take various specific forms, depending on the way 
dissimilarity is conceptualized. A precise example, from [13], is 
discussed in Section 4. In order for our results to remain as general 
as possible, we simply consider in this paper the class of C 2 
functions, i.e. whose first and second order derivative both exist 
and are continuous. 

And we assume they admit Hessian matrices that are nowhere 
negative semi-definite, i.e there is no admissible {P,}f =1 such 
that the Hessian of W ^{Pj}f =l ^ is negative semi-definite at 

{P,}f =1 . Weitzman's expected diversity function belongs to this 
class. It encompasses - but is not limited to - functions W with a 
positive definite Hessian matrix, i.e. that are strictly convex 
functions. 

Now let us take a step away from this initial metaphor, 
towards reality. Two modifications are brought into the 
formalism. First, rather than controlling directly the probability 
of survival P, of each species i=l,...,k, Noah can exert a 
protection effort within an admissible range, x, e [0,x,-], which 
is interpreted as the controlled increase of survival probability 
P, - say that x, is the increase of survival probability for species i 
resulting from a protection effort, e.g. an investment in a 
vaccination campaign, the provision of supplementary food, the 
protection and enhancement of habitat [6] . It is important to 
distinguish the effort from the change in the survival probability 
because P, is also determined by other factors, for there are 
ecological interactions among species. And this is where our 
second, most important, qualification appears: probabilities of 
survival are interdependent and the nature of those interactions 
are known. Nowadays, Noah can rely on the knowledge gained 
from the new and booming conservation biology literature on 
species distribution models and population viability analysis. See 
for instance [3], [14], [7], or [4] for a recent overview. Note that 
this literature does not take into account directly of species 
interactions; it just provides estimates of probabilities in space 
and time. From there, although applied econometric problems 
will have to be overcome, correlations between probabilities 
could be estimated. 

A group of experts can measure the marginal impact, say r,y,, 
that an increase in the probability of survival of a species h can 



have on the probability of survival of another species i. The 
experts can also appraise the impact of protection efforts on these 
probabilities. Assume, then, that the relationships between 
extinction risks are linear. Put differendy, a tractable approxima- 
tion of all those pieces of information can be summarized by the 
system (2) of linear equations: 

Pi = qt + x t + ^2 r 'h Ph, qi e [0,1 [ , x t e [0,x,-] . (2) 
h*i 

There are biological and economic factors that determines 
eligible efforts. Formally, admissible ranges of efforts are 
x f_ j [0,x,-] . Implicitly, additional efforts beyond the threshold X; 
have no effect on the survival probabilities. And we assume: 

Pi e Ui = [n,-,n,] C [0,1], Vf , V Xi e [0,3c,] . 

We denote P,- as the survival probability of species i without any 
conservation efforts, X; = 0,P,- > II, . In the absence of natural 
interactions, which corresponds to the case studied by Weitzman, 
we have rg, = 0,Vi,VA. A consequence is that in the very particular 
case with no ecological interactions and no conservation efforts, 
species i has a probability of survival q t . The survival probabilities 
interval, without ecological interactions, would thus take values 
ranging from P, = ^, to P, = P, +x,. 

Noah also has to cope with a budget constraint: 

X>,**,<P. (3) 
i=i 

where B is the total budget to be allocated to conservation - 
metaphorically, the size of the Ark - and /?,• is the cost per unit of 
effort to preserve species (. 

It is worthwhile making three remarks about this budget 
constraint. Firstly, it is assumed that changes in extinction 
probability are a linear function of expenditure. This may be 
inconsistent in real world applications where the marginal 
expense needed to reduce extinction risks is increasing. For 
example, [10] documents that the marginal preservation cost of 
threatened Australian birds increases when probability of 
extinction approaches zero. Weitzman rightly defends this 
linearity assumption as an acceptable approximation when the 
variation of probability falls in a sufficiently narrow range. But 
clearly, if costs are non linear and convex functions of efforts, an 
important qualitative result of our paper could change 
(Theorem 1 below may not hold any longer). Secondly, as a 
formal matter one could retrieve Weitzman's model with a 
simple change of variable, /?,■ = C,/AP, where C; is the cost per 
unit of increase of survival probability in the range 
AP, = n ; — n,. Thirdly, except when ecological interactions 
are negligible, Noah can increase the probability of survival of 
any species i via two different channels: a direct one by 
increasing the protection effort x,, at a cost yS, * X,- , and an 
indirect one through ecological interactions, due to the 
protection of another species j, with a cost j8 y - * X,-. 

Noah's Ark problem, when ecological interactions are taken 
into account, is then: 
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0,x,- 



w^jQ + u^p.jQ, (4) 



WoV(X) = W(V(X)), 



subject to (2) and (3). 

It will be convenient subsequently to work with matrix or vector 
expressions, written in bold characters. For any matrix M, let M T 
denote its transpose. Further, \ k is the (k x k) identity matrix, is 
the k dimensional column vector whose elements are all 1, and we 
recall the following definition of inequality between two k- 
dimensional vectors m and n with components m,- and ft, 
respectively: m^n if m,<w,, for all i= l,...,k. The other basic 
relationships between vectors are: f) m = n if »?, =«,-, for all 
i= \,...,k, ii) m < n if »?,<«,-, for all i = \,...,k, Hi) m<n if m,<n,-, 
for all i= \,...,k, and m^n. We also need basic matrix operations, 
"+", "-" and "*", that refer to, respectively the addition, the 
subtraction and the multiplication. 

Let us define: 



UoV(X) = U(V(X)). 



Under Assumption (INV), to each vector X corresponds a 
unique vector P = 'P(X). Therefore we can define Noah's problem 
with ecological interactions, the constrained maximization of a 
function of protection efforts X: 



max WoViX)- 
x 



UoV{X), 



subject to: 



ft T *X<B, 



(8) 



(9) 









■ 0 


rn • 


■ r xk ~ 




'Pi' 




"/V 




ft! 


, R = 




0 . 


■ r 2k 




Pi 


,P = 


h 




Ak. 




Jk\ 




• o . 




Pk. 




-h- 





'Pi' 




Pi- 




'xi ' 




'x\ ' 




Pi 




ll 








Xl 










x= 




X = 






Pk. 




-Pjc- 




x k _ 




Xk. 



In matrix form, the system (2) reads as: 

P = Q + X+R*P. (5) 

Throughout this article, we will assume: 
Assumption 1 (INV) The matrix \ k — R is invertible. 
Under Assumption (INV), the system (5) can be solved to give: 

P = A*(Q + X), (6) 

where A= [l*-R] _1 . 

Let V(X) = A * (Q + X) refer to the affine mapping from efforts 
to probabilities. Survival probabilities without protection policies 
are therefore: 

P = V(0*i k ), (7) 

where 0 * ( k is a vector made of k zeroes. Without ecological 
interactions, A is the identity matrix, P = Q and 
p = p + X = Q + X. 

Now we can plug (6) into (4) to get rid of probabilities, and 
express Noah's problem only in terms of efforts. Define the two 
composite functions, which here are mappings from the values 
taken by function V(X) to the set of real numbers: 



0*; k ^X^X. (10) 

Results 

Two questions arise: i) could anything general be said about the 
solution to the problem expressed by (8), (9), (10)? And if), taking a 
more practical stance, could we engineer a simple rule that 
approximates the general solution? 

Noah's policy is extreme 

Weitzman [13] showed that the solution to Noah's problem lies 
on the boundary of the efforts set. As the set of constraints is made 
of linear constraints, the boundary involves corners, e.g. x, = 0 or 
Xj = Xi, and possibly a segment between two corners, therefore 
with Xi e [0,Xi] for at most one species. This can be defined as an 
extreme policy. In words, the optimal protection policy gives full 
protection to a subset of species, partial protection for at most one 
species, and exposes the remaining species to the risk of no 
protection. 

But what if probabilities are interdependent? We show that 
when species interact, the optimal solution is also extreme. 

Theorem 1 The solution to Noah's Ark problem with ecological 
interactions, defined by (8), (9) and (10), is an extreme policy. 

Proof. The proof rests on two pieces of information: 

i) Noahs' problem is to maximize a continuous function over a 
compact set, therefore by Weiestrass extreme value theorem 
there exists a solution. 

ii) The Hessian matrix of W°V(X) + UoV(X) is not negative 
semi-definite, a statement we shall prove below. 

Item ii) violates the necessary second order condition for interior 
solutions to Noah's problem and, in combination with item i), 
leads to conclude the existence of a solution on the boundary of 
the efforts set. 

In order to prove item ii), because U °V(X) is linear, we just 
have to ensure that the Hessian matrix of WoV(X) is not negative 
semi-definite. Recall that V(X) is a ^-dimensional vector with 
typical element Th(X),h=\,...,k, and let Jp(X) stand for the 
Jacobian matrix: 
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Let us denote: 



/ P (X) = 



8Pi (X) 


8V\ (X) 


8Pi(X) 


dx\ 


5x 2 


8x k 


8V 2 (X) 


8V 2 (X) 


8V 2 (X) 


dx\ 


8x 2 


8xk 


8V k (X) 




8P k (X) 


dx\ 


8x 2 


8x k 



Note that, since each function T'a(X) is linear, the Jacobian 
matrix is made of invariant numbers, so we need not mention the 
application point X and we can simply refer to the matrix J-p. 

Denote V 2 WoV(X) the Hessian matrix of W°V(X), a k* k 
matrix with typical elements 8 2 ( W °V) / (dxtdxj) . From meticu- 
lous derivations of the composite function W°V(X), and after 
simplifications allowed by the linearity of the mapping V(X), one 
obtains: 



If V 2 W°P(X) is negative semi-definite, then for any nonzero 
vector m e R* we must have: 



m T * V 2 W°P(X) * m = m T * j£ * V 2 W(P) * J-p * m < 0 . 



Notice that J-p * m is simply a nonzero (k * 1) vector, which we 
may simply call n. Hence we can rewrite the above inequality as: 

n T *V 2 W(P)*n<0 , 

which would mean that V 2 H^(X) is negative semi-definite, a 
possibility that has been ruled out by assumption. ■ 



D, 



8W 



, u, 



p=p 

and define the two matrices: 



Di + Ui 
D2+U2 



8U_ 
~dP, 



A = 



D k +U k 



, Y = A * A. 



From simple calculations, the linearized problem in matrix form 
turns out to be: 



max Y * X + constant terms, 

x 



(11) 



subject to (9) and (10). 

As can be observed in the above approximation of Noah's 
problem, the introduction of ecological interactions changes the 
"slope" of the objective function to be maximized, which is now 
Y = A T * A instead of just A T . The crux, from the point of view of 
the present note, is to transform the information about ecological 
interactions conveyed by matrix R, into operational data via the 

matrix A=[l fc — R]~'. Given that I* — R is invertible, the 
computation of the matrix A is easily made and if A,y denotes a 
typical element of A, then Y is a ^-dimensional line vector of the 
type: 



Y=[ai , a 2 , ac k ] 



where 



A ranking rule for interacting species 

Theorem 1 is a qualitative result, that does not indicate which 
species should be granted protection and why. This brings us to 
our second question; it would be welcome to have an explicit and 
easy-to-use approximation of the general solution. Facing the same 
problem, this is the practical point of view adopted by [13], which 
he describes as "the main theme" of his paper (p. 1294). His 
formula (1) offers a ranking that is not really a solution to the 
original problem, but rather a first order approximation of an 
optimal policy. In order to achieve this, he replaces the objective 
function by its linear approximation. He then obtains a classical 
linear programming problem, whose solution is to assign grades Ri 
given by formula (1) to species (those grades depend on the model 
parameters) and order them in decreasing order of importance up 
to the point where the budget is exhausted. Those grades are 
exactly the practical ranking Noah is looking for. 

We follow the same approach here, i.e. we linearize the 
objective function. The astute reader knows that, in general, such 
approximations can be seriously misleading [2] and should not be 
followed blindly. Nevertheless, as proven in Theorem 2 below, 
there is something special about Noah's problem that makes this 
practice appropriate here. 



k 

a,= Y^{D h +U h )\ hi . 
h=l 

We can now define the "benefit"-cost ratios R' = a,-/ /?,- , or with 
explicit reference to relevant information: 

_. AD k 

tfEs—iJ^Dj+I^Aa, i=\,...,k- (12) 
L ' h=i 

As it is well-known, the argmax to the linear programming 
problem (11) is to fully protect the species with the highest grade 

R' , then the species with the second highest grade, and so on and 
so forth, up to the point where the budget is exhausted. It means 
that there exists a threshold value R* such that all species i with 

R' < R* are not embarked in the Ark, whereas those with grade 
larger than the threshold are all fully protected, except for at most 
one species with grade exactly equal the cutoff value R that is 
only partially protected. Let us call X w this policy, which can be 
described formally as follows: 
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if R' <R* , species i is granted zero protection, = 0 , 
if R' >R* , species i is granted full protection, =>x, = x, , 

< if R 1 = R* , species i is granted full protection (13) 
if there is enough budget, =>x; = X; , 

otherwise the remaining budget =>x; s ]0,x,[ . 
funds its partial protection, 



As shown in Theorem 2 below, X is a first order 
approximation of the optimal solution to Noah's Ark problem 
with ecological interactions. Put differently, there is a sense in 
which expression (12) can be taken for the new practical formula 
sought to construct in situ conservation priorities. Observe that the 
number assigned to each species / does not depend merely on its 
own "benefits" but actually on overall "benefits" generated by 

species i on all the species, ^\ _ ] (D/, + Uh) A;,,-, via ecological 
interactions. Therefore, a species with a strong own interest can be 
overridden by another, endowed with a less direct interest, but 
whose importance is enhanced because of its ecological role. Of 
course, when there are no ecological interactions, A is the identity 
matrix, with A„ = I, Ah, - = 0,V/i# i, and (12) boils down to Weitz- 
man's original system of grades for species i : 

—■ AP 

One can ask to what extent can we rely on formula (12) to build 
a hierarchy among species? Can a conservation policy be based on 
such an approximation? Baumol and Bushnell in [2] have 
famously attracted the attention on a number of potential flaws 
with linear approximations, two of them being important for the 
problem at hand: i) a linear approximation to a nonlinear program 
need not provide an answer better than a randomly chosen 
admissible answer, it) only if the objective function behaves 
monotonically in every variable within the admissible region can 
we be assured that a linear approximation will yield results which 
represent an improvement over the point where the linearization is 
made. Clearly, Noah's objective function does not meet this last 
condition, for an increase of the effort x, can improve the chances 
of species i at the expense of another species j (obviously so when i 
is a predator for j). 

Still, we can prove the following Theorem which establishes a 
special interest to the use of a linear approximation in this decision 
problem: 

Theorem 2 Consider the Noah's Ark Problem with ecological 
interactions, defined by (8), (9) and (10), and call X* its optimal 
solution. Then, 

i) the approximation of X* by X w , indicated in (13), offers an 
improvement compared to the absence of protection, 

ii) the approximation error, WoV(X*) — WoV(X ws ) , is no larger 

than K * (xV) 2 , where K= max{8 2 (W°V)/ (8 Xi dxj) }. 

Proof. Item i). The solution proposed in Theorem 2 is 
inspired from gradient methods used to find optimal solutions 



based on the property of iterative improvements, like the famous 
Frank-Wolfe algorithm. 

A first step is to replace the objective function by its first order 
Taylor approximation Z(X) computed at an admissible vector X° 
(here at the zero protection vector X° = 0 * i k ). Let us note VW°V 
the Gradient, a k * 1 vector with typical elements 8( WoV) / dxj , 
which corresponds actually to the vector Y = A T * A given in the 
text. 

Using those notations: 

Z(X) ~ WoV(0 * i k ) +V 1V°P(0 *i k ) T *(X-0*i k ). 

A second step is to find X that maximizes Z(X) subject to the 
relevant constraints. Since in Z(X) only the term 

VW°V(0 * i k ) X varies, this step is equivalent to maximize (1 1) 
subject to (9) and (10). And the policy X w presented in the 
Theorem 2 is exactly the maximizer of this linear programming 
problem. 

By definition of X w , we must have: 

Z(X w )>Z(0*i k ) . 

o VWoP(0*i k ) T *(X w -0*i k ) 
>S7WoP(0 * i k ) J *(0 * i k -0 * i k ) =0, 

so the vector X w — 0 * i k is an ascent direction for W°V. Although 
this means that the approximation Z(X) is non decreasing along 
this direction, it is not guaranteed that the non linear objective will 
behave similarly, i.e. we cannot yet conclude W°V(X W ) > 
WoV(0*i k ). 
By convexity of function WoV we can write: 

WoV (X w ) - WoV{0 *i k )>W WoV (0* i k ) T *(X w -0* i k ), 
and since we have established in (14): 



VW°V(0 * i k ) T *(X w -0* i k ) >0, 
we are led to conclude: 



WoV(X w )-WoV(0*i k )>0. 

Item ii) . Recall that V 2 WoV stands for the Hessian matrix of 
WoV. Using Taylor expansions, one can write: 

W°V(X*) = WoV(0 * l k ) 

+ VWoV(0 * i k ) T *X* + ^{X*Y *S/ 2 WoV{Z*)*X* , 
for some admissible vector Z*, and 
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WoV(X w ) = WoV(0 * i k ) +VWoV(Q * i k ) 

*X lv + 1 (X^) T *V 2 W°V(Z W ) * x w , 
for some admissible vector Z"'. Therefore 



WoP(X t )-WoP(X w ) =vw°p(o * i k ) T *(X* -x w 
+ ^ ! (X*) T *V 2 WVP(Z*) *x* 



-(x w Y*v 2 WoV{7r)*x v 



But, by definition of X w 



VW°V(0*i k ) *(X*-X w )<0, 



Pi- 



q2 + r 2 \q\ +x 2 + r 2 ixi 
\-r 2 \n 2 



(16) 



The grades also can be easily computed. They are: 



7? 1 - APl 



Di + Ux + r n (D 2 +U 2 )' 



R- 



AP 2 



\-r n r 2X \-r 2X r X2 



r u (Di + Ui) D 2 +U 2 



\-r\ 2 r 2 \ \-r 2 \r\ 2 



To further simplify, imagine that C\ = C 2 = C, APj = 
AP2 = AP. If ecological interactions are erroneously ignored, 
formally Noah assigns zero values by mistake to the system of 
interactions: r\ 2 = r 2 \ =0. Suppose, without loss of generality, that 
on this erroneous basis the first species ranks higher: 

R l = (D l + Ui)>R 2 = {D 2 +U 2 ) . 



W 



•V(X*) - W°V(X W ) < - (X*) T *V 2 W°V{T?) * X 



1 -(X W )\V 2 W°P{Z W )*X W 



<|[(x*)W] 2 -f [(x-)V] 



<k(x j *i k 



where K = max{ d 2 ( W°P) / (dxidxj) } . ■ 

The upper bound K for the approximation error mentioned in 
the above theorem if of course related to the non-linearity of 
W°V, formally captured by the second order derivatives 
d 2 (W°V) I '(dxidxj) . As a matter of interpretation, we can say 
that the stronger the curvature of W (the stronger preference for 
diversity if W is convex) the larger this upper bound. 



In other words D\ + U\=k* {D 2 + U 2 ), for some k> 1. 

Two questions arise. Could this ranking be reversed once 
interactions are properly taken into account? And, if the answer is 
affirmative, why? 

When the ranking is reversed: 



rUr 2 . 



D1 + U1 | r 2l (D 2 +U 2 ) ^ rnjJh + Ui) | D 2 +U 2 



\-r X2 r 2 \ \-r 2X r u 



\-r l2 r 2l \-r 2X r X2 



Since 1 — r\ 2 r 2 \ >0, and using D\ + U\ =k * (D 2 + U 2 ), the last 
inequality is equivalent to: 

<s> k + r 2i <kr X2 + \ , 



A Two-Species Example: Illustration and Discussion 

We close this note with an illustration using a simple two-species 
example. Let us first study to which extent the consideration of 
ecological interactions can alter priorities. Assume for simplicity 
that |ri 2 | < l,|f2i|< 1- The system (2) becomes: 



Pi 




11 




X\ 




p? 




+ 




+ 




CI?. 




x? 





0 r n 
rn 0 



k< 



1 



(since |ri 2 |< 1). 



So, a ranking reversal occurs when: 
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(17) 



Here the matrix I k — R is invertible since r \ 2 r 2 \ ^ 1. 
Solving the system of interactions: 



Pl = 



q\ +r\ 2 g 2 + x\+rux 2 
\-r\ 2 r 2 \ 



(15) 



In order to fix ideas, consider that k is arbitrarily close to one, 
i.e. the two species provide similar "benefits" and therefore a 
ranking reversal, if any, is due to the consideration of ecological 
interactions. Then note that for the above inequality to hold, 
necessarily r\ 2 > r 2 \ , which may occur in various interesting 
ecological configurations: 

i) Predation: species 1, a predator, feeds on species 2, its prey. 
So ^21 <0 whereas ri2>0. Giving conservation priority to 
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iii) 



the prey is the most effective way to enjoy the benefits of both 
species. 

Mutualism: for example plant-pollinator interactions, 
^12/21 >0. The synergistic relation between those two species 
is best enhanced by promoting species 2, which has the 
largest collective marginal impact. 

Competition: two species have to share a common resource in 
the same living area that cannot fully support both 
populations, hence r\2,f2\ <0, so conservation efforts focus 
on species 2 because its marginal negative impact is lower. 



Let us now examine the robustness of our results by specifying 
an expected diversity function. Denote G the number of genes 
joindy owned by the two species whereas Mj is the total number of 
genes owned by species i. Assume, as in [13] (expression (5)) that 
the expected (genetic) diversity function takes the following 
functional form: 



Of course, this is only an upper bound. In some cases, the 
approximation could also give the exact solution. To illustrate this, 
assume as before that C\ = C 2 = C, AP\ = AP 2 = AP, that utilities 
are identical, U\ = U 2 = U = 0, and the upper bounds on efforts 
are the same for the two species, X\ = x 2 = x. Assume also that the 
total budget can cover the protection cost of only one species, 
B = x * Cj AP. Noah then has to choose among two extreme 
policies, the first one (x\ = 0;x 2 =x) that provides the following 
expected diversity: 



W(0,x) = M l 



q\+r\2q2 + r\2X 

\-Y\2f2\ 



M- 



g2+?"21gl +X 

\-r\2fi\ 



M\ +rnq2 + rnx){q2 + r2\q\ +x) 



(\-rnr21) 

and the second one (x\ = x; x 2 = 0) with expected diversity: 



W(Pi ,P 2 ) = PiP 2 (M l +M 2 -G) + Pi(l- P 2 )M l 
+ P 2 (1-P 1 )M 2 + (1-P 1 )(1-P 2 )0 
= M l P l +M 2 P 2 -GP l P 2 . 

Considering relations (15) and (16) between efforts and 
probabilities, we obtain: 

WoP{x 1 ,x 2 ) = M ] P l {x u x 2 ) 

+ M 2 "P 2 (xi ,x 2 ) — GT\ (xi ,x 2 )V 2 {x\ ,x 2 ) . 



Two questions arise. Can we compare the true solution and the 
approximate solution? And can we estimate the error due to the 
approximation of the optimal solution? From Theorem 2, the 
upper bound on the error due to the approximation can be 
computed from the Hessian V WoV(X). In this two-species 
example, it is easy to derive the following formulae: 



d 2 ( WoV) /(dx l dx 2 ) = 8 2 ( W°V) /(8x 2 dxi) = - 



(l+ri 2 r 2 i) 
(l-r l2 r 2 i) 2 



G , 



d 2 (WoV)/(dx 1 y 



d 2 (WoV)/(dx 2 ) 2 



2r 2 \ 
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So the upper bound K on the approximation error, indicated in 
Theorem 2, is: 



W(x,0) = Mi gl + ? ' 2g2+ * + M 2 « 2 + r2 ^+ r2 ^ 
\-r n r 2 \ \-r\2r2\ 



-G 



(gi+n2g2+x)(ff2 + r2lgi+r2l*) 

(\-r\2r2\f 



It is optimal to protect species 2 if: 

W(0,x)>W(x,Q) 
o 

( 1 - mr 2 i )[M 2 ( 1 - r 2 \ ) - M\ ( 1 - r n )] > 
G[(l-r2i)(qi-q 2 + ruq2-r2iqi) + (n2-r2i)x] 



In the particular case where G = 0, then M t = 8W /dPi = Dj, 
and the above condition boils down to a very simple expression: 

W(0,x) > W(x,0)o %± = ^ < , 
M 2 P> 2 i—r\2 

a condition which is also necessary for the approximated solution 
to select species 2 (remember condition (17)). It comes as no 
surprise that the optimal solution and its approximation concide, 
since when G = 0 the upper bound on the approximation error is 
zero, as can be seen from expression (18). 
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K = {\-r X2 r 2 i)- 2 *G*max{-2r 21 ,-(l+n2r2i),-2r 12 } , (18) 

a value which depends only on the number of genes owned jointly 
by the two species, G, and on the ecological interaction terms, rjj. 
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